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ABSTRACT 

Experimental  results  of  synthesis  of  some  real  textures 
using  a  class  of  spatial  interaction  models  known  as  simul¬ 
taneous  autoregressive  (SAR)  models  are  given.  The  role  of 
statistical  inference  methods  in  texture  modeling  is  empha¬ 
sized. 
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1.  Introduction 


Texture  is  an  important  characteristic  for  the  analysis 
of  many  images.  Image  texture  can  be  described  in  terms  of 
a  two-level  mechanism,  one  to  describe  the  basic  primitives 
composing  the  texture  and  the  other  to  describe  the  spatial 
layout  of  the  primitives.  The  primitives  may  be  of  deter¬ 
ministic  shape  such  as  squares,  circles,  hexagons,  or  may  be 
of  completely  random  shape.  Textures  such  as  ceiling  tiles 
and  bricks  can  be  grouped  under  the  first  category  while  tex¬ 
tures  such  as  grass  and  sand  fall  under  the  second  category. 

It  may  be  possible  to  model  patterns  such  as  bricks  and  tiles 
by  using  the  placement  rule  model  [Rosenfeld,  Troy,  1970; 
Zucker,  1974] .  Such  a  description  may  be  very  complex  for 
random  textures  such  as  cork  and  grass  due  to  the  random  shapes 
and  orientations  of  the  primitives. 

We  take  the  view  that  the  given  texture  is  generated  by 
an  underlying  random  field  characterized  by  appropriate  para¬ 
metric  models.  One  of  the  prime  characteristics  of  this 
approach  is  that  the  models  considered  for  the  textures  are 
generative,  i.e.,  given  the  structure  of  the  model  and  reason¬ 
able  estimates  of  the  parameters  characterizing  the  model  syn¬ 
thetic  textures  close  to  the  original  textures  may  be  obtained. 
Several  generative  stochastic  models  are  known  in  the  litera¬ 
ture  for  texture  analysis.  By  appending  the  successive  rows. 
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autoregressive  and  moving  average  models  (ARMA)  and  their 
variants  {McCormick,  Jayaramamurthy,  1974] ,  and  by  seasonal 
autoregressive  models  [Delp,  et  al.,  1979].  Since  we  are 
dealing  with  two-dimensional  textures,  intuitively  it  would 
be  preferable  to  represent  the  textures  using  2-0  models. 
Two-dimensional  causal  ARMA  models  have  been  suggested  for 
textures  in  [Tou,  1980].  However,  the  causal  models  are 
only  a  restricted  set  of  models. 

Recently,  a  different  class  of  spatial  interaction  models 
known  as  conditional  Markov  (CM)  models  have  been  used  for 
modeling  textures  [Hassner,  Sklansky,  1978A;  1978B;  1980; 

Cross,  1980].  Hassner  and  Sklansky  [1978A;  1978B]  have  pointed 
out  the  equivalence  between  the  Gibbs  field  and  llarkov  random 
field  and  generated  some  synthetic  binary  textures.  The 
coding  method  [Besag,  1974]  for  obtaining  the  estimates  of 
the  CM  models  has  been  briefly  mentioned.  Cross  [1280]  has 
generated  a  wide  variety  of  binary  and  eight-level  synthetic 
textures  using  the  Monte  Carlo  simulation  procedure  [Metro¬ 
polis,  Rosenblutri,  et  al.,  1953].  The  Monte  Carlo  procedure 
is  iterative  in  nature  and  involves  considerable  computational 
time  for  generating  synthetic  textures  with  large  numbers  of 
gray  levels. 

Our  main  emphasis  in  this  paper  is  to  illustrate  the  appro¬ 
priateness  of  SAR  models  for  textures.  Although  under  Gaussian 
situations,  there  exists  a  correspondence  between  the  classes 


of  SAR  and  CM  models  [Besag,  1974]  in  that  given  a  SAR  model, 
an  equivalent  CM  model  can  be  found,  this  correspondence 
is  true  only  in  second-order  properties.  Hence  for  non- 
Gaussian  situations,  the  two  models  are  different.  In  the 
equivalent  SAR  and  CM  models,  usually  the  SAR  models  are 
characterized  by  a  lesser  number  of  parameters.  The  main 
attraction  of  the  SAR  models  considered  here  is  that  the 
methods  can  be  easily  generalized  to  non-Gaussian  situations. 
Such  a  generalization  to  non-Gaussian  situations  using  CM 
models  is  unattractive  due  to  the  computationally  expensive 
iterative  schemes  used  for  synthesis  [Cross,  1980]. 

The  organization  of  this  paper  is  as  follows:  In  Section 
2,  we  give  the  theoretical  variograms  of  several  SAR  models. 
Many  of  these  variograms  exhibit  periodic  patterns  found  in 
the  variograms  of  natural  textures.  Section  3  discusses  the 
role  of  statistical  inference  methods  in  texture  modeling 
using  SAR  models.  Specifically,  we  review  the  role  of  spe¬ 
cific  estimation  methods  and  model  selection  rules  discussed 


in  [Chellappa,  1980  ;  Chellappa,  Kashyap,  1981B;  Chellappa, 
1981].  In  Section  4,  we  give  the  experimental  results  of 


fitting  SAR  models  to  cork,  grass,  wood,  sand  and  paper  taken 


from  Brodatz's  album  [Brodatz,  1956].  The  quality  of  fit  is 
evaluated  using  visual  inspection  and  empirical  variograms 


of  the  reconstructed  textures, 
given  in  Section  5. 


Finally,  a  brief  discussion  is 
F _ a  ■  r - - - 
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2.  Appropriateness  of  SAR  models  for  texture  synthesis 
We  intend  to  illustrate  the  usefulness  of  SAR  models 
for  textures  in  two  steps,  using  the  synthetic  patterns 
generated  by  known  SAR  models  and  the  theoretical  vario- 
grams  of  SAR  models.  In  an  earlier  report  tChellappa, 

1930] ,  several  synthetic  patterns  were  generated  using 
SAR  models.  Contrary  to  the  existing  belief  [Modestino, 
et  al. ,  1979]  that  spatial  interaction  models  such  as  SAR 
models  are  incapable  of  exhibiting  local  replication 
attributes,  an  essential  ingredient  of  texture,  several 
of  the  synthetic  patterns  generated  by  the  SAR  models  do 
show  local  replication. 

Another  way  of  judging  the  appropriateness  of  SAR 
models  is  by  using  the  variogram.  Suppose  {y(s),s€ft}, 
ft-{s-(i, j) ,  lsi,j<M}  are  the  observations  from  the  texture. 
Then  the  variogram  at  displacement  (k,2.)  is  defined  as 

V(k,»)  -EKy(s)  -  y(s+(k,t))I2l 

-  <2-1' 

where  ^  is  the  normalized  autocorrelation  function.  An 
expression  for  the  theoretical  variogram  may  be  obtained  by 
assuming  that  the  given  set  of  observations  {y(s)>  obeys  a 
finite  lattice  SAR  model  in  (2. 2) -(2. 3)  [Kashyap,  1980A; 

1S30B;  1981] .  Prior  to  that  we  partition  the  finite  lattice 
ft  into  two  sets  ft^  and  ftfi,  defined  below: 


and 


=  (s=(i,j):  s€G  and  (s+r)€Ji> 


°I  -  0  -  °B 


The  representation  of  {y(s)}  is  different  in  and  as 
given  below: 

y(s)  *  I  0  y  (s+y)  +  /fa)(s),s€£i  (2.2) 

y€N  T 


and 


y  (s)  =  E  8  y.(s+y)  +  /Iu)(s),s€Q  (2.3) 

Y€N  Y  1  B 


Wiiere 

y1(s+(k,i))  with  s=(i,j) 

=  y(s+(k,£))  if  (s+(k,£)Hft 
*  y  ( (i+k-l)mod  M+l,  (j+£-l)mod  rt+1]  if 

(s+(k, *))/«. 

In  (2.2)  and  (2.3),  {u>(.)}  is  an  independent  and  identically 
(IIu)  distributed  noise  sequence.  The  set  N  is  known  as  the 
neigiibor  set,  and  depending  upon  the  specific  choice  of  N, 
causal,  semicausal  and  noncausal  models  are  obtained  [Jain, 
Jain,  1978] .  Causal  or  more  generally  unilateral  neighbor 
sets  are  finite  subsets  of  the  half  plane  S+  defined  in 
[Goocunan,  Ekstrom,  1980].  The  observation  set  {y(s)}  is  not 
Markov  with  respect  to  N,  i.e., 


p  (ys)  |  all  y  (y) ,  Y5*s) 

7*  P (y  (s)  |  all  y(s+Y)>  Y*N)  (2.4) 


as  is  true  for  CM  models 


Due  to  the  specific  finite  lattice  representation, 


the  normalized  autocorrelation  function  pk  ^  at  lag(k,Jt) 
can  be  written  as  [Chellappa,  Kashyap,  1981A] 


E  exp  /=I  X  n[(s-l)T(k,Jt)]/|  |  ys  |  |2 
s€fl  u 

E  1/t |Ps|  j2 
s€fl 


(2.5) 


where 

y  s  =  (i-eT*8)#8€fl,x0-^-  (2.6) 

ip  =  Col.  [exp(/^T(s-l)Ty)  ,y€N]  (2.7) 

S 

and 

T 

1  =  (1,1) 

The  variograms  V(k,&)  related  to  pk  z  as  in  (2.1)  are  plotted 
in  Figure  1  for  several  SAR  models.  The  details  of  the  raoaels 
are  given  in  Table  1.  For  each  model  we  have  plotted  vario¬ 
grams  along  tne  four  directions  (0,i) , (i,0) , (i,i) , (i,-i)  in 
a  discrete  lattice.  The  structure  of  the  possible  variograms 
with  two-dimensional  SAR  models  is  quite  varied  and  many  of 
them  possess  oscillatory  behavior,  a  characteristic  of  tne 
variograms  of  natural  textures  with  periodic  patterns  [Schacnter, 
et  al. ,  1978] . 
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I 
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mating  the  parameters  of  the  model  given  the  structure 
of  the  model.  For  SAR  models  with  unilateral  neighbor  sets 
the  classical  least  square  (LS)  estimates  are  consistent 
and  efficient  for  the  Gaussian  case  but  the  LS  estimates 


are  not  consistent  for  nonunilateral  neighbor  sets  [Ord, 
1975;  Kashyap,  Chellappa,  1981] .  To  obtain  consistent  and 
efficient  estimates,  maximum  likelihood  (ML)  estimates  can 
be  obtained  by  making  appropriate  assumptions  regarding  the 
distribution  of  { co  ( s )  >  -  Due  to  the  Jacobian  of  the  trans¬ 
formation  matrix  from  (u>(.)}  to  { y  ( - )  }  not  being  unity  for 
non-unilateral  neighbor  set  SAR  models,  the  log  likelihood 
function  for  Gaussian  {y(s)}  is  non-quadratic  in  the  para¬ 
meters.  This  requires  the  use  of  computationally  expensive 
gradient  schemes  like  Newton- Raphson  to  obtain  ML  estimates. 
An  iterative  scheme  with  reduced  computations  has  been  deve¬ 
loped  in  [Chellappa,  1980;  Chellappa,  Kashyap,  1981B; 
Cnellappa,  1981]  which  yields  approximate  ML  estimates.  The 
approximate  ML  estimates  W,]}  are  obtained  as  limits  of 


6t+l'  ®t  9iven  below: 


9t+i  -  <?  -  ^  ?)“1(V  -  j-  U),t=0,l,2,3 . 


(2.8) 


(2.9) 


3t  *  m"2  n(y(s)  ”  ®t5(s))2'  t=o,i,2,3,... 


where 

S  =  E  z(s)zT(.s),  mxm  matrix  (2.10) 

a  ~ 

U  =  I  z(s)y(s),  m*l  vector  (2.11) 

a  " 


V  *  E  C  ,  mxl  vector  (2.12) 

~  a  ~s 

R=E(SS  -  C  CT) ,  mxm  matrix  (2.13) 

a  ~s“s  ~s~s 

Cg  =  Col.  [cos  ^-(s-l)T  y,y€N] 

and 

Sg  =  Col. [sin  ^-(s-l)Ty,y€N] 

The  initial  value  0Q  is  chosen  as  follows: 

0„  =  S_1U. 

~  U  ~  ~ 

All  the  summations  in  (2.8)- (2.13)  are  over  sea  and 
m  is  the  dimension  of  0.  The  use  of  approximate  ML 
estimates  in  texture  modeling  has  been  emphasized  in  [Chellappa, 
1980;  Chellappa,  Kashyap,  1981B]  using  synthetic  patterns. 

The  quality  of  the  regenerated  synthetic  textures  using 
approximate  ML  estimates  is  better  than  the  synthetic  tex¬ 
tures  using  inconsistent  LS  estimates. 

Thus  far  in  our  discussions  we  have  assumed  that  the 
underlying  structure  N  of  the  model  is  known.  In  practice, 
the  appropriate  N  is  to  be  estimated  from  the  given  texture 
data.  Asymptotically  consistent,  transitive  and  parsimonious 
decision  rules  to  choose  the  appropriate  N  for  the  given 


I 


i 


texture  have  oeen  given  in  [Chellappa,  1980;  Kashyap, 
Chellappa,-  19£>1]  .  The  derivation  of  these  decision  rules 
may  be  found  in  [Kashyap,  et  al.,  1981].  The  relevance 
of  these  decision  rules  in  texture  modeling  has  been  illus¬ 
trated  in  [Chellappa,  1980;  Kashyap,  Chellappa,  1981]  using 
synthetic  patterns.  Several  SAR  models  were  fitted  to  the 
data  generated  oy  a  known  SAR  model  and  the  test  statistic 
gn  given  in  (2.14)  was  computed  for  each  model: 

g  =  M2logf  -  Zlog(l-20jc 
n  n  ~n~sn 

+  finin'  +  ^  (2-14) 


where 


9sn 


(  £ 
s€Q 


T 

C  C 
~sn~sn 


T 

+  S  S  ) 
~sn~sn' 


(2.15) 


The  model  corresponding  to  the  minimum  of  gn  is  chosen  as 
the  appropriate  model.  This  decision  rule  correctly  eliminates 
several  inappropriate  models.  One  way  of  judging  the  in¬ 
appropriateness  of  a  model  is  by  evaluating  the  quality  of 
the  synthetic  texture  generated  by  the  model.  The  synthetic 
textures  corresponding  to  clearly  inappropriate  models  are 
not  similar  to  the  original.  But  the  synthetic  patterns  corres¬ 
ponding  to  SAR  models  whose  neighbor  sets  include  the  neighbor 
set  of  the  true  model  are  very  similar  to  the  original  making 
visual  judgment  subjective.  Due  to  its  preference  for  parsi¬ 
monious  models,  the  decision  rule  minimizing  gn  correctly 
eliminates  the  over-parameterized  models. 


4.  Real  texture  synthesis  using  SAR  models 

Our  discussions  so  far  have  been  concerned  with  syn¬ 
thetic  image  patterns  generated  by  known  SAR  models.  The 
next  step  is  to  investigate  the  appropriateness  of  SAR 
and  CM  models  for  natural  textures  like  sand,  cork,  grass, 
ecc.  For  our  experiments,  we  chose  5  textures  from  Brodatz's 
album  [1956] .  The  64x64  textures  selected  were  cork,  sand, 
grass,  wood  and  paper.  We  report  the  results  of  fitting 
SAR  models  to  the  textures  mentioned  above  as  a  sequence  of 
experiments . 

In  the  first  experiment,  different  SAR  models  were  fitted 

to  cork  using  the  estimation  scheme  in  (2.8)— (2.9) .  The 

details  of  the  fitted  models  are  in  Table  2;  the  estimate  6 

and  the  test  statistic  g  in  (2.14)  are  in  Table  3.  The 

n 

patterns  corresponding  to  the  fitted  models,  generated  using 
a  Gaussian  pseudo-random  number  generator,  are  given  in 
Figure  2.  The  synthetic  textures  (1,2),  (1,3),  and  (2,1) 
corresponding  to  SAR  models  Ngi»NS2'  an^  ^gg  are  not  900<i*  *n 
fact,  models  Ng^  and  Ng2  give  non-stationary  patterns.  The 
quality  of  the  synthetic  textures  (2, 2) -(3, 2)  is  significantly 
improved  when  neighbors  farther  than  the  nearest  are  included, 
as  in  models  Ng4  through  Ng7  in  Table  2.  Since  the  reduction 
in  variance  is  not  significant  by  increasing  the  number  of 
parameters  from  10  to  14,  the  decision  statistic  gn  for  the 
10  parameter  model  Ng4  is  less  than  those  for  the  12  and  14 
parameter  models  Ngg  and  Ng7. 


However,  the  18  parameter  model 


N^7  is  preferred  to  the  lower  order  models.  Observe  that 
the  patterns  corresponding  to  the  last  four  models  are 
very  similar. 

The  synthetic  cork  textures  corresponding  to  models 
NS4  ~  NS7  di-sPlay  the  diagonal  pattern  present  in  the  ori¬ 
ginal  cork.  However,  the  synthetic  textures  seem  to  have 
more  abrupt  variations  in  their  intensities  than  the  original 
corx,  i.e.,  distinct  dark  and  white  patches  are  present  in 
cue  synthetic  textures  rather  than  smooth  variations  as 
in  the  original  cork. 

wote  that  the  synthetic  textures  in  Figure  2  were  gener¬ 
ated  using  a  Gaussian  pseudo-random  number  generator.  Better 
quality  can  be  achieved  if  the  noise  sequence  used  to  generate 
the  synthetic  textures  is  derived  by  using  the  histogram  of 
the  actual  residuals.  The  actual  residuals  (uj(.)}  were  pro¬ 
duced  by  using 

«  =  |[H(0)y  -  al]  (2.16) 

wnere  B,9"  and  a  are  the  estimates  of  6,9,  and  a,  respectively. 
We  used  three  standard  programs  BDC0U1,  MDGC  and  GGVCR  from 
the  IMSL  package  in  that  order  to  generate  a  set  of  random 
deviates  whose  histogram  approximates  the  histogram  of  the 
exact  residuals.  The  exact  residuals  were  tallied  into  twenty 
bins  using  BDC0U1.  The  normalized  frequencies  obtained  from 

BDC0U1  were  fed  as  inputs  to  MDGC.  The  cumulative  distribution 
function  of  the  residuals  obtained  from  MDGC  was  used  to 


generate  the  required  random  deviates.  The  results  of  syn¬ 
thetic  generation  using  the  histogram  matched  random  numbers 
are  given  in  Figure  3.  Figure  (1,1)  is  the  original  cork 
texture,  and  window.  (1,2)  is  the  synthetic  texture  using 
the  14  parameter  SAR  model  and  Gaussian  random  number  genera¬ 
tor.  The  synthetic  texture  using  the  histogram  matched  resi¬ 
duals  is  in  window  (1,3) . 

The  results  of  fitting  SAR  models  to  the  other  textures  under 
consideration  are  given  in  Figure  4.  The  model  fitted  was 
14  parameter  SAR  model  Nsg  and  histogram  matched  residuals 
were  used  to  generate  the  synthetic  textures.  The  synthetic 
paper  texture  (1,2)  in  Figure  4  retains  the  vertical  direc¬ 
tionality  present  in  the  original  in  (1,1).  As  mentioned  be¬ 
fore,  some  large  patches  not  in  the  original  texture  are  pre¬ 
sent  in  the  synthetic  image.  The  horizontal  streak-like  pat¬ 
tern  present  in  the  original  wood  texture  (2,1)  is  reasonably 
picked  up  in  the  synthetic  wood  (2,2),  with  a  few  additional 
patches.  The  synthetic  grass  (2,3)  corresponding  to  the 
natural  grass  in  (1,3)  is  not  as  good  as  cork  or  wood.  This 
is  probably  due  to  the  fact  that  grass  appears  to  be  more 
inhomogenous  compared  to  cork  or  wood.  As  in  the  case  of  paper 
and  cork  the  directionality  present  in  sand  (3,1)  is  picked  up 
in  the  synthetic  texture  (3,2).  Another  way  of  judging  the 
quality  of  texture  synthesis  is  to  see  if  the  second  order 
properties  such  as  the  variogram,  etc.  of  the  synthetic  texture 


are  close  to  the  original  data.  To  answer  this  query,  vario- 
graias  of  original  textures  and  synthetic  textures  are  given 
side  oy  side  in  Figure  5.  The  synthetic  textures  correspond 
to  SAxi  model  Ngg  and  were  generated  using  histogram  matched 
resiuuals.  Note  that  in  all  the  cases  the  variograms  of  syn¬ 
thetic  textures  are  reasonably  similar  to  the  original  vario¬ 
grams.  One  that  is  worthy  of  pointing  out  is  the  variogram 
of  the  synthetic  wood  along  (0,i). 

It  is  probable  that  a  64x64  data  set  may  be  large  for 
a  14  parameter  model.  Large  images  may  be  blocked  into  small 
images  and  each  block  may  be  modeled  separately  [Delp,  et  al., 
1979].  The  results  of  fitting  12  parameter  symmetric  SAJR 
models  to  32x32  blocks  of  four  textures,  cork,  sand,  paper, 
and  grass,  are  given  in  Figure  6.  The  synthetic  textures  were 
constructed  using  histogram  matched  residuals.  Compared  to 
the  64x64  synthetic  cork  texture,  the  32x32  synthetic  cork  in 
window  (2,2)  has  fewer  abrupt  intensity  variations.  The  sand 
and  paper  textures  also  show  such  results.  However,  the 
32x32  synthetic  grass  (1,4)  is  not  as  good  as  the  64x64  syn¬ 
thetic  grass  in  Figure  6. 


5.  Discussion 


We  have  given  some  experimental  results  of  our  attempts 
to  synthesize  a  class  of  real  textures.  The  textures  con¬ 
sidered  belong  to  the  class  of  microtextures.  There  have 
been  several  earlier  attempts  to  synthesize  real  textures, 
notably  the  work  reported  by  Cross  [1980]  and  Garber  and 
Sawchuck  [1981] .  Cross  considered  a  wider  class  of  textures 
than  attempted  in  our  study.  In  the  Monte  Carlo  procedures 
used  for  synthesis  purposes,  the  size  of  the  state  set  in¬ 
creases  rapidly  with  the  number  of  gray  levels.  Hence,  to 
avoid  excessive  computations,  two  and  eight  level  images 
obtained  by  an  equal  probability  quantization  procedure 
[Haralick,  et  al.,  1973]  have  been  fitted  with  CM  models. 

Coding  estimates  have  been  used  to  generate  synthetic  tex¬ 
tures.  Good  results  have  been  reported  for  the  textures 
considered  in  our  study.  It  has  been  concluded  that  CM 
models  are  poor  fits  to  textures  such  as  brick,  water,  and 
wood.  It  is  not  clear  if  the  bad  fits  are  due  to  inefficient 
coding  estimates  or  due  to  the  distinctly  inhomogeneous  struc¬ 
ture  of  wood,  water  and  brick  due  to  the  reduction  in  number 
of  gray  levels  or  due  to  the  inadequacy  of  the  CM  models  them¬ 
selves.  Our  limited  experiments  with  SAR  models  have  given 
reasonably  good  results  for  wood  and  it  appears  that  much  larger 
parameter  sets  are  required  for  water.  Due  to  the  extremely 
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high  regularity  present  in  the  brick,  probably  a  two-stage 
procedure  may  be  appropriate  [Cross,  1980,  p.  137],  In  the 
first  stage  we  can  model  the  line  patterns  for  bricks  and  sub¬ 
sequently  fill  in  the  areas  using  a  SAR  model. 

A  number  of  synthesis  techniques  have  been  suggested  in 
[Garber,  Sawchuck,  1981]  for  binary  and  gray  level  textures. 
Mostly  unilateral  models  have  been  considered  compared  to 
the  more  general  SAR  models  used  in  our  study.  Consequently, 
LS  estimates  have  been  used  to  generate  synthetic  images. 

No  quantitative  rules  are  available  in  their  study  for  choos¬ 
ing  an  appropriate  model.  Good  results  have  been  reported 
for  binary  images  and  some  gray  level  images.  Better  results 
might  have  been  obtained  on  smaller  images  since  the  number 
of  parameters  required  would  be  much  smaller  leading  to  better 
accuracy  in  estimation. 
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Table  1.  Details  of  SAR  Models  Whose  Variograms 
are  in  Figure  1.  0  is  assumed  to  1.1111. 


Model 

Identification 
Fig.  la 

Fig.  lb 


Fig.  lc 


N 

Parameter  0 

(-1,0), (-1,1),  (0,1),  (1,1) 

0_1  0*.18,9_^  -.  1806 

0Q  1=1.1011,91  -1.0390 

(-1,0), (-1,1), (0,1), (1,1) 

e-i,oss0i,o“*5256 

(1,0) , (1,-1) , (0,-1) , (-1,- 

U  0_i  1=91  _1=  -.2480 

VrV-r*5081  , 

9i,i"e-i,-i"  ~*2876 

(-1,0) , (0,-1) , (-1,-1) , 

9-i,o=i-0388' 

(-2,0) , (0,-2) 

V-r-9046 

0_1  _1=-.7288 

0-2,l=S  "'18O4,0O,-2= 

-.1048 

(-1,0), (-1,1), (0,  1) 

9-l,0=9l,0®’28' 

*-* 

O 

** 

1 

H* 

O 

1 

H* 

9o,-r0o,iss--14 

9-l,l=0l,-l=*22 

Fig.  Id 


Model 

Identification 


Fig.  le 


Table  1 
Ccont. ) 


N 


Parameter  9 


(-1,0) , (0,-1) , (-1,-1) 


9-i,o“-9704' 

90,-l“*9735' 
9_1  _x=  -.9686 


e-l,l“0l,-l"#23, 

9o,r9o,-r  ‘*14 


Fig.  If 


(-1,1) , (0,1) , (1,1) 
(1,-1) , (0,-1) , (-1,-1) 


91,1  9-l,-l 
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models  fitted  to  cork 


Model 

Number 

Symmetric  Neighbor  Set 

NS1 

(0,1) , (1,0) 

NS2 

(0,1) , (1,0) , (1,-1) , (-1,-1) 

NS3 

(0,1), (1,0) , (1,-1) , (-1,-1) , (0,2), (2,0) 

NS4 

(0,1), (1,0), (1,-1) , (1,1) , (0,2) , (2,0) 

(-2,1) , (2,1) , (1,2) , (-1,2) 

NS5 

Ns4  U  {(2,2),  (-2,2)} 

WS6 

Ns5  U  { (3,0) ,  (0,3) } 

*^S7 

Ns6  U  {(1,3), (3,1),  (-1,3), (-3,1)} 

Table  3 

Estimate  |3  anti  test  statistics  g  for  SAR  models 
fitted  tc  cork.  Variance  of  Sata  =  714.04 


Model 

Number 

r 

8 

NS1  ! 

200.02 

Non-stationary 

NS2 

169.78 

Non-stationary 

NS3 

168.52 

23423.0 

NS4 

157.73 

23225.0 

NS5 

157.68 

23243.0 

NS6 

157.71 

23257.0 

NS7 

155.28 

23207.0 

The  theoretical  variograms  of  the  SAR 
models  in  Table  1.  V(0,k)  versus  k 

is  plotted  in  the  first  row,  V(k,0) 
versus  k  in  the  second  row  ,  V(k,k) 
versus  k  in  the  third  row  and  V(k,-k) 
versus  k  in  the  fourth  row. 


Figure  2. 


Synthesis  of  cork  texture  using  SAR  models. 
(1,1)  is  the  original  cork,  (1 , 2 ) , (1 , 3) , (2 , 1) 
were  generated  by  models  N_,  -  N§3  of  Table  2. 
Patterns  (2 , 2) ,  (2 , 3) ,  (3 , 1) “and  (3,2)  were 
generated  by  Ng4  -  Ng7. 


Figure  3.  Results  of  using  histogram  matched  residuals. 

(1,1)  is  original  cork,  (1,2)  was  generated  by 
Ngg  of  Table  2  using  Gaussian  pseudorandom 
numbers  and  (1,3)  was  generated  by  N  ,  and 
histogram  matched  residuals.  Sb 


Figure  4.  Results  of  fitting  SAR  models  to  other  textures. 

All  the  synthetic  textures  were  generated  using 
Ngg  and  histogram  matched  residuals. 

(1.1)  and  (1,2):  original_and  synt.^cic  paper. 

(1,3)  and  (2,3):  original  and  synthetic  grass. 

(2.1)  and  (2,2):  original  and  synthetic  wood. 

(3.1)  and  (3,2):  original  and  synthetic  sand. 


Original  Cork 


Figure  5.  Variograms  of  natural  and  synthetic 
textures.  Synthetic  textures  were 
generated  using  SAR  model  Ng6  and 
histogram  matched  residuals. 


Figure  6.  Results  of  fitting  SAR  models  to 

32x32  blocks  of  textures.  All  the  synthetic 
textures  were  generated  using  Ng5  and 
histogram  matched  residuals.- (1,  Jj  and 

(1.2) :  original  and  synthetic  paper. 

(1.3)  and  (1,4):  original  and  synthetic 
grass.  (2,1)  and  (2,2):  original  and 
synthetic  cork.  (2,3)  and  (2,4):  original 
and  synthetic  sand. 
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